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Introduction 


The NOAA/NASA Pathfinder program was created by the Earth Observing System 
(EOS) Program Office to determine how existing satellite-based data sets can be 
processed and used to study global change. The data sets are designed to be long 
time-series data processed with stable calibration and community consensus 
algorithms to better assist the research community. 

The Ocean Altimeter Pathfinder Project involves the reprocessing of all altimeter 
observations with a consistent set of improved algorithms, based on the results from 
TOPEX/POSEIDON (T/P), into easy-to-use data sets for the oceanographic 
community for climate research. This report describes the processing schemes used 
to produce a consistent data set and two of the products derived from these data. 
Other reports have been produced that: a) describe the validation of these data sets 
against tide gauge measurements and b) evaluate the statistical properties of the 
data that are relevant to climate change. 

The use of satellite altimetry for earth observations was proposed in the early 1 960s 
(‘Williamstown Report,” MIT, 1970). The first successful space based radar altimeter 
experiment was flown on SkyLab in 1974 (McGoogan, et al., 1974). The first 
successful satellite radar altimeter was flown aboard the Geos-3 spacecraft between 
1975 and 1978 (The Geos-3 Project, 1979). While a useful data set was collected 
from this mission for geophysical studies, the noise in the radar measurement and 
incomplete global coverage precluded it from inclusion in the Ocean Altimeter 
Pathfinder program. This program initiated its analysis with the Seasat mission, 
which was the first satellite radar altimeter flown for oceanography (Seasat, 1 982). 
This mission was followed by the U.S. Navy’s Geosat satellite in 1 985 to 1 989 
(Geosat , 1987, 1990a, 1990b). The European Space Agency’s ERS-1 spacecraft 
was flown from 1991 to 1996 (ERS-1 System, 1992). ERS-1 was followed by the 
joint U.S. and French T/P mission from 1992 to the present (TOPEX/Poseidon, 1994, 
1 995) and ESA’s ERS-2 spacecraft in 1996 (ERS-2 Spacecraft, 1995). Table 1 
summarizes the dates of these missions, the data coverage, and the length of time 
for which a Pathfinder data set has been constructed. 

The second section of this report describes the algorithms that are used to process 
these data sets into a consistent cross mission data set for climate studies. The third 
and fourth sections describe two data products that have been produced from the 
corrected data. In section 3, a collinear product is described that is, essentially, time 
series of sea surface height along the track of the spacecraft. In these data only 
minimal interpolation and smoothing have been done to the measurement. In section. 
4, a grid data set is described where monthly anomalies relative to an annual 
reference surface have been created from ail missions to provide a decade long data 
set of sea level anomalies. 

The validation of the data sets are described in Pathfinder Report #2 and a statistical 
description of sea surface topography based on these data are described in Report 
#3. 
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• Table 1: Satellite altimeter missions and data sets between 1975 and 2000. 

Mission Dates Geophysical Data Pathfinder Data 

Record Dates Dates 


Seasat 

6/26/78-10/10/78 



1 7 day near repeat 


7/7/78-9/10/78 

7/78-9/78 

3 day repeat 


9/11/78-10/10/78 

9//78- 10/78 

Geosat 

3/12/85-12/31/89 



Primary Mission 


3/30/85 - 9/30/86 

4/85-9/86 

Exact Repeat 
Mission 


11/8/86-12/31/89 

11/8/86 - 8/22/89 

ERS-1 

7/17/91 -2/6/96 



Phase A - 3 day 

commissioning 

phase 

7/31/91-12/20/91 

7/31/91 - 12/20/91 

Not processed 

Phase B - 3 day 

12/28/91 -3/30/92 

12/28/91 -3/30/92 

12/28/91 -3/30/92 

Phase C - 35 day 

4/14/92-12/20/93 

4/14/92-12/20/93 

4/14/92-12/20/93 

Phase D - 3 day 

12/23/93 - 4/10/94 

12/23/93-4/10/94 

12/23/93-4/10/94 

Phase E - 1 68 day 

4/10/94-9/27/94 

4/10/94 - 9/27/94 

4/10/94 - 9/27/94 

Phase F - 168 day 

9/27/94-3/21/95 

9/27/94-3/21/95 

9/27/94-3/21/95 

Phase G- 35 day 

3/21/95-6/2/96 

3/21/95-6/2/96 

3/21/95-6/2/96 

TOPEX/POSEIDON 

8/10/92 - Present 

9/23/92 - Present 

9/23/92 -Present 

ERS-2 - 35 Day 

4/20/95 - Present 

4/21/95 - Present 

4/21/95 - Present 























The satellite altimeter measurement of surface height is observed by merging two 
complicated measurements, the radial position of the space-based ranging 
instrument relative to the center of the planet, and the round-trip travel time of a 
radar or light pulse sent from the spacecraft (see Figure 1). To transform satellite 
altimeter range and orbit measurements into accurate sea-surface elevation data, 
a variety of models, algorithms, and corrections need to be adopted and applied. 
These include the determination of the satellite position within a consistent 
geodetic reference frame, the correction for atmospheric range refraction and 
radar ranging corrections, and the removal of unwanted geophysical effects, such 
as tides and atmospheric loading. Most of the adopted PATHFINDER 
corrections and algorithms are consistent across satellite missions, but some are 
specific to an individual mission. Table 2 provides a summary of these many 
parameters for the satellites listed in Table 1 . In this section, these algorithms are 
described. 

• Figure 1 A schematic of the satellite altimeter measurement of sea surface 
topography. 
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Table 2 Summary of Algorithms used to process satellite altimeter data. 


II 1 li 1 lii L, Ml 

mJ.'* 

TOPEX 

POSEIDON 

ERS-1 

Geosat 

Seasat 

II.A. Position 






Orbit 

T/PPOD2 

JGM-3 

T/P POD 2 
JGM-3 

Delft 

DGM-E04 

Goddard 

JGM-3 

Goddard 

JGM-3 

Geodetic ref 

ITRF 

ITRF 

ITRF 

ITRF 

ITRF 


T/P 

T/P 

T/P 

T/P 

T/P 







II. B. 

Instrument 






Sea State Bias 

Walsh 

Caspar BM4 

Gaspar BM1 

Gaspar BM3 

Gaspar 

Retracking 

None 

ONES 

ESTeC 

Goddard 

Goddard 

Calibration 

H 

ONES 

Range bias 

Cal mode and 
oscillator drifts, 
Range bias 

Range bias 







II.C. 

Atmospheric 

Refraction 






Dry Atmos. 

Saastamoinen 
w/ ECMWF 

Saastamoinen 
w/ ECMWF 

Saastamoinen 
w / ECMWF 

Saastamoinen 
w/ ECMWF 

Saastamoinen 
w / FNOC 

Wet Atmos. 

Measured 

Measured 

Measured 

nggsBHHj 

Measured 

Ionosphere 

Measured 

DORIS 

IRI95 

IRI95 

IRI95 







il.B. 

Geophysical 

Models 






Earth Tide 

Modified Wahr 
(1981) 

Modified Wahr 
(1981) 

Modified Wahr 
(1981) 

Modified Wahr 
(1981) 

Modified Wahr 
(1981) 

Ocean Tide 

Schrama & Ray 

(1994) 

SR960104 

Schrama and 
Ray (1994) 
SR960104 

Schrama and 
Ray (1994) 
SR960104 

Schrama and 
Ray (1994) 
SR960104 

Schrama and 
Ray (1994) 
SR9601 04 

Load Tide 

SR9601 04 

SR960104 

SR960104 

SR960104 

SR960104 

Pole Tide 

IERS 

IERS 

IERS 

IERS 

IERS 

Geoid 

Gradient 

UTCSR 
Mean Surface 

UTCSR 
Mean Surface 

UTCSR 
Mean Surface 

UTCSR 
Mean Surface 

UTCSR 
Mean Surface 

Atmos Load 

IB w/ECMWF 

IB w/ECMWF 

IB w/ECMWF 

IB w/ECMWF 

IB w/FNOC 
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2.1. Orbits 


2.1.1. TOPEX^POSEIDON Precise Ephemerides 

The precise satellite ephemerides for T/P are identical to those used in the 
second generation T/P Geophysical Data Records (GDRs) produced by the 
TOPEX Project Office. These ephemerides are of unprecedented accuracy, and 
they are, in large part, responsible for the enhanced accuracy of T/P over any 
other historical or present-day satellite-altimeter mission. The ephemerides are 
produced by the Space Geodesy Branch at NASA Goddard Space Flight Center 
(GSFC). They are computed for nominal 1 0-day arc segments, corresponding to 
the ground-track repeat period of the satellite. They are based on a combination 
of laser and doppler tracking data and on force models derived from the JGM-3 
earth gravity model, a T/P-based ocean-tide model, and a 'box-wing' satellite 
model for drag and radiation forces. The first-generation of these orbits are 
described by Nerem et al. (1993) and Tapley et al. (1 994). The present version of 
the orbits is described by Marshall et al. (1 995), who examine in detail the error 
budget and present evidence for a radial orbit accuracy of about 2 to 3 cm rms. 

The Joint Gravity Model 3 (JGM-3) is an update to the JGM-1 and JGM-2 models 
developed at GSFC and the University of Texas and described by Nerem et al. 

(1 994). The JGM-1 model was developed before the launch of T/P. It was the 
result of a multi-year effort to improve the earth's gravity model by a new inversion 
of tracking data on over 30 satellites, altimeter data from Seasat and Geosat, and 
direct gravity measurements on the earth's surface (land and marine gravimetry). 
The JGM-2 model was a 'tuning' of JGM-1 after the launch of T/P by inclusion of 
1 50 days of T/P tracking observations. The JGM-2 model was used for 
computing the first-generation orbits on the original T/P GDRs. The JGM-3 
model represents a further tuning of JGM-1 by inclusion of more tracking data on 
T/P, and especially the inclusion of about 40 days of GPS tracking. The JGM-3 
model is described in detail by Tapley et al. (1996). 

The tidal perturbations on the satellite are deduced from a new ocean-tide model 
developed by Ray et al. (1 994) that was produced primarily from T/P altimetry. 
Poleward of T/P's maximum 66 p latitude, the model was supplemented by 
Schwiderski's older hydrodynamic model. Spherical harmonics through degree 
and order 15 are used from this ocean model to compute satellite perturbations, 
except for the resonant terms (degree-2, order-2 semidiumals, degree-2, order-1 
diumals) which are taken from the tracking-based GEM-T3S model of Lerch et al. 
(1992). In addition to the ocean tide, of course, the earth tide also produces 
satellite perturbations. The earth model developed by Wahr (1 981 ) is used. Any 
error in the earth tide will be absorbed, and hence adequately compensated for, 
by the GEM-T3 ocean resonance terms. 

The box-wing model which is used for modeling forces arising from atmospheric 
drag, solar and earth radiation, and spacecraft radiation has been described in 
detail by Marshall and Luthcke (1 994a; 1994b). 

Four different tracking systems are available on T/P: a) French Doppler 
Orbitography and Radiopositioning Integrated by Satellite (DORIS) system; b) 
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Satellite Laser Ranging (SLR); c) satellite-satellite tracking via the Tracking and 
Data Relay Satellite System (TDRSS), and d) Global Positioning System (GPS). 
For these computed orbits, only the two 'operational' tracking methods (DORIS 
and SLR) are used. (The TDRSS and GPS data have been used by GSFC and 
other groups for orbital accuracy assessments and for other studies.) During the 
force-model integration and the fitting to the tracking data, one-per-revolution 
acceleration adjustments are also included to handle any remaining modeling 
errors; the theory for this adjustment is described by Colombo (1 986) and Cretaux 
et al. (1994). Many other details concerning these orbits, including critical details 
concerning the geodetic reference frame and tracking station positioning, are 
described in the references mentioned above. 

Analysis of collinearsea surface height residuals reveals dominant mesoscafe 
oceanographic features as well as providing a statistical assessment of orbit 
models, geophysical and environmental height corrections. No adjustments 
were performed to correct for systematic orbit errors for the maps shown in 
figures 2-4. 

• Figure 2 T/P Sea Surface Height Variability 1993-1995. 
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2.1.2. ERS-1 Precise Ephemerides Version 2.0 

The precise satellite ephemerides for the ERS-1 mission (Phases B to G) have 
been computed at the Delft Institute for Earth-oriented Space Research (DEOS), 
Delft University of Technology. The computations used the DGM-E04 gravity 
field model (Scharroo et al., 1 997). The DGM-E04 is a tuning of the JGM-3 
model (Tapley et al., 1 996) especially for handling ERS-1 and ERS-2. The 
tracking data on ERS-1 consists of satellite laser ranging only. 

The radial orbit accuracy for the Phase C orbits (eighteen 35-day repeat cycles) is 
approximately 6-8 cm rms. The average ERS-1 altimeter crossover residual rms 
is about 1 0 cm for 1 -sigma, including an estimate for geographically correlated 
orbit errors due to gravity. 

• Figure 3 ERS-1 Phase C Sea Surface Height Variability. 
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2.1.3. Geosat Precise Ephemerides 

The satellite ephemerides for Geosat were computed at GSFC (Williamson and 
Nerem, 1 994). They are based on TRANET Doppler tracking data collected at 
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some 45 ground stations and on force models that include the JGM-3 Earth 
gravity model and a T/P-based ocean-tide model. 

The tidal forces on the satellite are computed using the same ocean-tide mode! 
that was used in the T/P ephemerides (Ray et al., 1 994). Spherical harmonics 
through degree and order 15 are used from this ocean model to compute satellite 
forces, except for the resonant terms (degree-2, order-2 semidiumals, degree-2, 
order-1 diumals) which are taken from the tracking-based GEM-T3S ocean-tide 
solution of Lerch et al. (1992). 

• Figure 4 Geosat Sea Surface Height Variability During the First Year of the 
Exact Repeat Mission. 
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2.1.4. Seasat Precise Ephemerides 



Precise Seasat ephemerides based upon JGM-3 were also computed at GSFC. 
These orbits are based upon TRANET Doppler, Unified S-Band (USB), and 
satellite laser ranging data. The force modeling is as was described above for 
Geosat. The USB tracking station coordinates required for Seasat were not 
available from the JGM-2/3 solution, so these were recovered from analyses of 
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the Seasat tracking data, holding fixed the TRANET sites and most of the lasers. 
Fifteen USB site locations were adjusted, as were the three SAO laser sites at 
Arequipa, Mt Hopkins, and Natal. These SAO site locations were adjusted 
because, with the current orbital accuracies, laser data biases in the early part of 
the Seasat Mission were both more detectable and separable from the site 
locations. 

The predominant relevance of the USB tracking is that, in the early part of July 
1 978, TRANET tracking is not available, and there is only very sparse laser 
tracking. Thus the recovery of USB coordinates enables the computation of 
these earlier orbits in the same reference frame as the latter orbits. 

it should be noted that the USB data do have some associated unmodeled 
ionosphere errors and so the earlier orbits are expected to be somewhat less 
accurate than in the latter period. 


Z2. Reference Frames 

2.2.1. Geodetic Reference Frame 

The geodetic reference frame, which is consistently used in the PATHFINDER 
reprocessing for all altimeter satellite data, is the International Terrestrial 
Reference Frame (ITRF). The ITRF definition and other current modeling 
recommendations of the International Earth Rotation Service are documented in 
their “IERS Conventions,” (McCarthy, 1996). 

In the context of reference systems and station positioning, it is interesting to note 
the special case of Geosat, which was tracked only by the TRANET Doppler 
system. Were it not for the Seasat mission, which employed both TRANET and 
SLR tracking, it would be very difficult to adjust Geosat TRANET tracking sites 
into the ITRF. Thus, Seasat is important not only for its altimetry, but also for its 
tracking data. Appropriate use of those data allows Geosat to be tied into the 
same geodetic reference system as T/P. Work is still in progress to strengthen 
this tie. A consistent geodetic reference frame is crucial to climate-related 
studies. 


2.2.2. Reference Ellipsoid 

The reference ellipsoid used for all Pathfinder altimeter products is consistent with 
the official T/P GDRs. 


Table 3 Reference Ellipsoid Parameters. 


Term 

Value 

Equatorial radius 

6378136.3 m 

Flattening 

1/298.257 

Potential 

62636858.702 rffVs 2 

G Me 

398600.4415 km 3 /s 2 ~ 
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2.3. instrument Corrections 


2.3.1. Sea State Bias 

The sea-state bias correction compensates for the bias of the aftimeter range 
measurement toward the troughs of ocean waves. This bias is thought to arise 
from three interrelated effects: a tracker bias, a skewness bias, and an 
electromagnetic (EM) bias. The correction for this bias is dependent on some 
percentage, usually some variable percentage, of the significant wave heights. 

2.3.1 .1. TOPEX 

For the NASA TOPEX altimeter we adopted the EM-bias model devised by Ed 
Walsh (T/P GDR Handbook, Callahan 1993) with no additional tracker bias. This 
model is parametrically equivalent to the one used on the original TOPEX GDRs, 
but the values themselves are different for the following reasons: 

1 . Backscatter coefficients (o 0 ) have been adjusted for the calibration results of 
Callahan, Hancock, and Hayne (1994; this reference gives adjustments only 
through cycle 40; other cycles are personal communications). 

2. Wind-speed algorithm implementation, see below. 

The present algorithm is thus 

Eg. 1: TIP SSB = Hsw (0.0029 + 0.0038 U - 0.00015 U 2 ) 

where, Hsw is the altimeter observed significant wave height in meters, U is the 
altimeter observed 10-m wind speed in meters/second inferred from the adjusted 
(j 0 by using the modified Chelton-Wentz (MCW) algorithm with an additional -0.63 
dB offset. Witter and Chelton (1991) describe the MCW algorithm. Callahan, 
Morris, and Hsiao (1994) describe the offset. The T/P SSB is in meters. 

2.3.1.2. POSEIDON 

The Poseidon altimeter is the CNES altimeter aboard the T/P satellite, nominally it 
is turned on for every 1 0 th cycle of T/P observations. For this altimeter we 
adopted an updated version of the algorithm published by Gaspar et al. (1994), 
'BM4'. (Update coefficients from P. Gaspar, personal communication.) This is a 
four-coefficient model, determined by an empirical study of Poseidon crossover 
residuals. The model uses as independent variables the observed significant 
wave height H sw and the inferred wind speed U, the latter being deduced from the 
backscatter cross-section coefficient, a 0 , by use of the MCW algorithm. The 
algorithm is 

.Eg. 2P0SEID0NSSB = H sw (0.047+ 0.0023 U- 0.000112 U 2 - 0.001 H sw ) 
Where, H sw is in meters, U in meters/second, and POSEIDON SSB is in meters. 
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2.3.1. 3.ERS-1 


For the ERS-1 altimeter we adopted a simple sea-state bias correction using a 
constant percentage of significant wave height. This percentage is taken as 
5.5%, following an unpublished study of ERS-1 crossover data by Philippe 
Gaspar. The fact that the percentage (5.5) is much higher than the nominal 2% 
used for some altimeters indicates a significant tracker bias, similar though not 
quite as large as that seen on Seasat. Gaspar found little statistically significant 
dependence on wind speed, although this dependence has been established for 
other altimeters and in ground-based measurements. 

In summary the correction used in this Pathfinder product is 

. Eq. 3 ERS-1 SSB = 0.055 H sw 
Where, H sw is in meters, and ERS-1 SSB is in meters. 

2.3.1.4.Geosat 

For correcting the Geosat altimeter for the sea-state bias, we have relied on the 
study by Gaspar and colleagues (1996, Pathfinder Report #4) who recommend: 

. Eq. 4 Geosat SSB = H sw (-0.025 - 0.00145 U + 0.00020 H sw U) 

Where, H sw is in meters, U is the inferred wind speed in meters/second, and 
Geosat SSB is in meters. The wind speed U is deduced from the backscatter 
coefficient (a 0 ) using the Freilich and Challenor (1994) Rayleigh-based relation. 
This gives the wind speed at 1 9.5 m. The speed is divided by 1 .057 to obtain a 
1 0-m neutral stability wind speed. 

2.3.1 .S.Seasat 

Prior estimates of the sea state bias for SEASAT ranged between 6-7% of the 
significant wave height (Bom et al. 1982, Douglas and Agreen, 1983). After 
retracking the ocean returns, our preliminary estimate of the sea state bias is 
2.7% of the significant wave height (P. Gaspar, personal communication, 1997). 
An empirical study of the SEASAT crossover residuals is currently being 
performed by Phillipe Gaspar to provide a formal estimate of the sea state bias. 

2.3.2. Range Retracking 

2.3.2. 1 .TOPEX - Carried out by Ernesto Rodridquez of JPL but we 
have not applied these to the data set 

2.3.2.2. POSE1DON - Carried out and applied to GDR by CNES 

2.3.2.3. ERS-1 - Carried out and applied to GDR by CERSAT 

2. 3. 2.4. Geosat and Seasat 


16 



•altimetry returns 


The Ice Altimetry Group of NASA/Goddard Space Flight Center (GSFC) and 
Hughes SIX have developed algorithms for netracking polar ice altimetry returns 
to account for the inability of the on-board tracker to follow the topography of the 
ice sheets. The Seasat and Geosat altimeter data returns from the sensor data 
records (SDR’s) were retracked employing the GSFC Version 2 retracker. 

Details of the retracking concepts and the Version 2 algorithm are described 
below and can also be accessed at: 

http://crevasse.stx.com/ia home/retrack/qsfcretrackdoc.960725.html 
GSFC Retracking Concept 

The GSFC algorithms are based on the concept that the return can be 
represented mathematically using a function described by Parsons 1979 over the 
Gaussian ocean surface. Typical ocean returns are shown in Figure 5. The 
measurement corresponding to the midpoint of the ramp gives the median 
elevation within the pulse-limited footprint. The on-boand tracker tries to align the 
return so the midpoint of the ramp occurs at the tracking gate. For the ocean 
surfaces this works very well as indicated by the retracking correction varying 
from 3-6 cm. The longer vertical line indicates the position of the on-board 
tracking gate, whereas the short vertical line indicates the midpoint of the ramp, 
where the tracker should have been tracking to, as determined by the retracking 
process. 

♦ Figure 5 Typical Ocean Waveforms over Antarctic Ocean from Geosat 



range gate 


An iterative Bayesian least squares procedure is used to fit the digitized waveform 
to the appropriate function. For the single ramp return, a modified version of the 
ocean return function is used where an additional slope parameter is added after 
the ramp. 
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The mid-point of the ramp, Beta(3), is chosen as the correct range point, and the 
retracking correction is the displacement of that mid-point from the tracker gate. 
The lineartrailing-edge of the function simulates the slowly decreasing return 
from a diffuse reflecting surface as the pulse expands over the beam-limited 
footprint. A nine-parameter function is used to simulate the double-peaked 
returns. 
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Before fitting, a priori values for the coefficients are calculated from a linearly- 
filtered return. A set of standard deviations is assigned to these a priori values 
indicating how well each coefficient could be approximated from the filtered 
return. The digitized waveform consisting of digitized "gates" are also weighted in 
the solution to try to put the heaviest emphasis on fitting to the ramp location (s) 
which defines the retracking correction. The function coefficients are then solved 
for using the unfiltered digitized return by iterating until the parameters defining 
the ramp positions, Beta(3) and Beta(6), converge to within two percent. After the 
3rd iteration some of the standard deviations on the parameters are decreased to 
stop the procedure from either diverging or going into a limit-cycle situation. This 
process usually works, but it is using a linear solution to a non-linear problem and 
sometimes unexpected results occur. When the functional fit process fails, the 
retracking correction is calculated using the threshold-type algorithm discussed 
above for near-specular returns. 
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2.3.3. Calibration Issues 


2.3.3.1 .TOPEX 

Calibration of the TOPEX satellite altimeter is an ongoing concern. This 
PATHFINDER project uses the most recent calibration adjustments. We have 
applied corrections based on internal calibration-mode estimates of the altimeter 
range drift, provided by the Wallops Flight Center PATHFINDER team members. 
We have also corrected for the error in the oscillator drift that was discovered in 
July 1996. 

In addition we have applied a correction for drift in the measured backscatter 
coefficient, o 0 , although this affects estimates of sea-surface height only through 
the EM-bias correction. The a 0 adjustment (up through T/P cycle 40) is described 
by Callahan, Hancock, and Hayne (1994). 

2.3.3.2. POSEIDON 

POSEIDON range measurements are corrected on-board the spacecraft with an 
internal calibration. As of now, there is no reason to suspect a drift in the 
Poseidon altimeter aboard the T/P satellite after the internal calibration has been 
applied by CNES. This statement rests primarily on comparisons of the Poseidon- 
based sea-surface elevations with elevations measured at oceanic tide gauges. 

2.3.3.3. ERS-1 

There is a known bias in the ERS-1 altimeter of between 40 and 45 cm. The 
altimeter appears to be ranging short, so derived sea-surface heights are too high 
by this amount. The bias has been estimated from observations at the ERS-1 
Venice calibration site (Francis et ai., 1 993) and from crossover comparisons with 
T/P altimetry (Le T raon et al., 1 995; Brian Beckley, unpublished results, 1 997). 
Our colinear dataset has NOT been corrected for this bias, but the grid data 
product is adjusted (see Table 7). 

2.3.3.4. Geosat 

There are several corrections that have been applied to the Pathfinder Geosat 
products to compensate for some long-term drift-like errors. Two of them, a 
calibration-mode adjustment to the altimeter and a drift adjustment to the 
oscillator, are shown in Figure 6. 

The 'cal-mode‘ (or calibration-mode) drift can be detected by the altimeter itself 
when a synthetic pulse is passed through part of the internal electronics (see 
MacArthur et al., 1987). The result of such processing is shown in the top portion 
of the figure. In addition to a very small long-term drift, there are anomalous 
'humps' that occur near the beginning of each calendar year and are correlated to 
on-board spacecraft temperatures. We have evidence that these humps are real, 
because they tend to "straighten out" estimates of global mean sea-level 
computed from the altimetry. Therefore, we have applied this correction to the 
Pathfinder products in the form shown. 
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Surface Height Additive Surface Height Additive 

Correction for Oscillator Drift, cm Correction from Cal Mode, cm 



The Geosat clock or oscillator was found to drift about 3 parts per billion per year. 
This was noticed during the original processing of the altimetry, and it was 
'corrected' midway through the mission by a shift of 10 parts per billion. We have 
unapplied this shift and then corrected for a constant drift of roughly 2 mm per 
year in sea level. The bottom portion of Figure 6 shows the correction that we 
have applied, including the allowance for the shift during October 1987. 


• Figure 6 Variation in Geosat clock oscillation. 


GEOSAT ERM Surface Elevation Corrections 



GEOSAT ERM Surface Elevation Corrections 



2.3.3.5.Seasat 

Height calibration values from Seasat are not reported herein. Upon completion 
of a formal estimate of the sea state bias by Gaspar, geocenter offsets relative to 
T/P will be computed and reported in Table 7 at 
http://neptune.asfc.nasa.gov/oceans.htmi . 
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2 A Atmospheric Range Refraction Corrections 


2 A1 . Dry-Atmosphere Range Refraction Correction 

The dry atmospheric range correction is designed to correct the observed 
altimeter sea-surface elevations for propagation delay caused by refraction 
through the dry component of the atmosphere. It requires an estimate of the 
barometric surface pressure at the satellite nadir point. Given this pressure P in 
mbar, the dry-atmosphere correction is approximated by (Saastamoinen, 1972) 

• Eq. 5 Dry Atmosphere Correction = -2.277*P*[1.0 + 0.0026*cos(q)] mm 


Where, q is twice the latitude. This approximation assumes that the dry 
atmosphere obeys the perfect gas iaw and that it is in hydrostatic equilibrium. 
The small second term shown in brackets accounts for the variation in surface 
gravity owing to the Earth's ellipsoidal shape. 

The pressure data used here are extracted from the operational model of the 
European Centre for Medium-Range Weather Forecasting, which is updated 
every 12 hours (except for T/P where a special 6-hour product has been used). 

• Figure 7 Distribution and size of the dry atmosphere range refraction 
correction (mm), during the ERS-1 mission Phase C cycle 2, i.e. during the 35 
day interval beginning on 19 May 1992. The mean correction during this 
period was 2305 mm, with a standard deviation of 26 mm. 



2210 2245 2280 2315 2350 
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2.4.2. Wet-T roposphere Range Refraction Correction 

The wet troposphere range correction is designed to correct the observed 
altimeter sea surface elevations for the propagation delay caused by refraction 
through atmospheric water vapor and cloud liquid water. Although it is smaller 
than the corresponding dry atmosphere range correction, the wet correction is 
more complex, with possibly rapid variations in both time and space. This 
correction can vary from just a few millimeters in dry, cold air to more than 40 cm 
in hot, wet air. 

T/P: The correction is based on the measurements of the T/P on-board 
microwave radiometer. This is a 3-channel (operating at 1 8, 21 , and 37 
Gigahertz), nadir-viewing instrument, which provides nearly direct measurements 
of the wet-troposphere range refraction correction by monitoring the strong water- 
vapor absorption line centered at 22.235 GHz. The radiometer and its 
corresponding data-processing algorithms are described in a series a three 
papers published in the IEEE TGRS (Ruf et al., 1 995; Janssen et al., 1 995; 

Keihm et al., 1 995). The radiometer was calibrated in flight during the first few 
months of the T/P mission; this calibration work is described by Ruf et al. (1 994). 
Further work concerning the validation of this instrument includes Stum (1994) 
and Keihm and Ruf (1 995). The accuracy of this calibration is now being called 
into question, in 1997, because of a suspected 1-2 mm per year drift in the height 
measurement that may be caused by the radiometer drift. 

ERS-1: The correction is based on the measurements of the on-board 
microwave radiometer. This is a 2-channel (operating at 23.8 and 36.5 
Gigahertz), nadir-viewing instrument. The radiometer and its corresponding data 
processing algorithms and initial performance are described by Bernard et al 
(1993). 

Geosat: Unlike some altimeter missions, this mission did not include an on- 
board microwave radiometer that would normally monitor the nadir atmospheric 
water content. For Geosat wet-troposphere corrections, we must rely on other 
sources. Earlier versions of Geosat datasets have used climatological analyses 
from Nimbus-7 SMMR data, operational products from the Fleet Numerical 
Oceanography Center, and data from the TIROS Operational Vertical Sounder 
(TOVS) and from the SSM/I instrument on the DMSP satellite. The last two are 
based on the work of Emery et al. (1 990). 

For this Pathfinder product, we use the computed 'precipitable water" field that 
has been produced in the recent Reanalysis Project of the NOAA National 
Centers for Environmental Prediction. For a description of this project, see 
Kalnay et al. (1 996). The columnar water vapor has been converted to a 
troposphere range refraction correction by use of the high-order algorithm of Liu 
et al. (1 990); this algorithm takes account of the variability of atmospheric 
equivalent temperature. 
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• Figure 8 Distribution and size of the wet atmosphere range refraction 
correction (mm), during cycle 30 of the T/P mission, i.e. during the 10 day 
interval starting 1 July 1993. The mean correction during this time was 
156mm, the maximum 454mm and the minimum Omm. 



2.4.3. ionosphere Range Refraction Correction 

The ionosphere range refraction correction to the altimeter data compensates for 
the propagation delay in the travel time caused by the presence of free electrons 
in the atmosphere, primarily the ionosphere. 

2.4.3.1.TOPEX 

For the first time, the TOPEX altimeter allows for a direct measurement of this 
effect, rather than relying on numerical models as in previous satellite-altimeter 
missions. This is done through the use of a dual-frequency altimeter, one 
frequency in the C band (5.3 GHz) and the other in the traditional Ku band (13.6 
GHz). Because the ionosphere refraction delay is, to a good approximation, 
inversely proportional to frequency squared, the two altimeter height 
measurements yield a direct first-order correction for the time delay. Further 
discussions of the TOPEX ionosphere measurements can be found in Imel 
(1994). 

The raw ionosphere correction in TOPEX measurements is slightly noisy and 
requires smoothing before it is applied to the altimeter range measurements. 

This is done by a local Gaussian convolution filter, with a width of 20 one-second 
altimeter measurements. Figure 9 shows a long section of altimetry (beginning 
approximately at MJD 49062.8) with the raw ionospheric corrections and their 
smoothed counterpart in red. The ionosphere's characteristic bimodal 
enhancement about the equatorial electrojet, its dependence on local time, and 
its dependence on solar declination are all apparent in Figure 10, which displays 
the mean ionosphere correction during the first two months of 1 994 at four 
different local times. As anticipated, the maximum correction occurs at mid-day. 
Also, since the sun is in the southern hemisphere during January and February, 
the ionosphere correction is enhanced in that hemisphere. 
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The ionospheric range delay computed from the TOPEX data has an expected 
accuracy of 0.5 cm. 


• Figure 9 Section of T/P altimetry (beginning at approximately MJD 49062.8) 
with the raw ionospheric corrections and their smoothed counterparts in red. 
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• Figure 10 The mean ionosphere correction (mm) during the first two months 
of 1994 at four different local times . Mean local times of measurements given 
above each plot. 



2. 4.3.2. POSEIDON 

For the Poseidon altimeter, this correction is based on ionospheric 
measurements made with the DORIS dual frequency receiver aboard the T/P 
satellite and DORIS beacons at many ground locations. As the DORIS receiver 
overflies the beacons, the Doppler shift difference in the two frequencies is 
monitored and then combined with the Bent ionospheric model (Bert and 
Lleweilyn, 1 973) to generate a map of the ionosphere corresponding to the 
appropriate local time. 
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When the TOPEX altimeter is operational, its dual-frequency ionospheric 
measurements may be compared to the estimated DORIS measurements. 

Figure 1 1 shows an example, based on one series of measurements beginning at 
approximately MJD 49062.8. The raw unsmoothed TOPEX measurements of 
ionospheric path correction are compared to the estimated correction for DORIS, 
shown in blue. 

• Figure 1 1 The raw unsmoothed TOPEX measurements of ionospheric range 
refraction are compared to the estimated correction from DORIS model , 
shown in blue. Beginning at approximately MJD 49062.8 



Time aioag track (seconds) 


2.4.3.3.Geosat and ERS: IRI95 Ionosphere Correction for Altimeters 

When lacking direct altimeter measurements of the ionospheric delay we rely on 
other data or on models of the ionosphere. For ERS-1 , Seasat, and Geosat 
altimetry, we use the International Reference Ionosphere 1 995 (IRI95), 
specifically Version 1 3 dated January 1 996, in combination with a special 
ionosphere-effective, solar index. This index is obtained from ionosonde 
measurements, which during the Geosat time period comprised about 70 globally 
distributed stations. 

The IRI95 model represents a marked improvement over earlier models (Bilitza, 

1 997, Bilitza et al., 1 995). However, by comparing it to actual ionosphere 
observations made by the TOPEX altimeter, it is easy to find periods when the 
IRI95 model yields an inaccurate correction for altimetry. Figure 1 2 is an 
example. Near the equatorial electrojet, the IRI95 underestimates the required 
correction by order 1 0 cm. (Also shown for comparison is the ionosphere 
correction from the DORIS system, as used for POSEIDON altimetry.) We 
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emphasize that this plot has been purposely selected as an example of inferior 
performance by the IRI95 model; it is not typical. 


• Figure 12 Comparison of TOPEX TEC measurement, DORIS ionosphere 
model and IRI95 models. 


Ionospheric Range Correction 



2.5. Geophysical Confections 

2.5.1 . Solid-Earth Body Tide Correction 

The Earth's body tide is modeled as a purely elastic response to the lunar and 
solar tidal potentials. The algorithm used here for all satellites is identical to that 
used for the T/P GDR data, in order to maintain consistency with the ocean-tide 
models, which by and large depend on T/P data. 

The lunar and solar tidal potentials are computed from the semi-numerical 
expansion of Cartwright and Tayler, described by Cartwright and Edden, (1973), 
extrapolated linearly to the 1 990 epoch. The complete expansion is used, with 
terms of both second and third degrees in the potential. Only the first term, from 
the permanent tide, has not been used. By agreement (within the T/P Science 
Working Group), the permanent tide is included in the geoid model. 

The body-tide vertical displacement depends on the adopted h 2 and h 3 Love 
numbers. For the T/P and Pathfinder projects, the following Love numbers were 
adopted: 

h 2 = 0.609 for all terms of degree 2, except 
= 0.52 for the Ki tide and its nodal sidelines. 
h 3 = 0.291 for all terms of degree 3. 
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These terms are approximately those computed by Wahr (1981 ), with the 
exception that Wahr found h 2 = 0.606 for long-period tides and a much more 
complex behavior in the diurnal band surrounding the nearly diurnal free-wobble 
frequency. Also, all latitude dependence in the Love numbers has been 
neglected. As noted above, however, we have followed the T/P algorithm for 
consistency with the ocean tide models 

2.5.2. Ocean Tide Correction Version 2.0 

The ocean tide removal in the altimeter surface height measurement is the largest 
of all standard corrections. For example, in an analysis of collinear differences of 
sea-surface heights, Ray, Koblinsky, and Beckley (1 991 ) found that the ocean 
tides were responsible for more than 80% of the signal variance. 

The Pathfinder ocean tide model for all satellites is actually a combination of 
many models. The primary deep-ocean model is an updated version of Schrama 
and Ray (1994). This is version 9601 04 (hereafter SR9601 04), but it is nearly 
identical to version 950308, which was studied extensively by Shum et al. (1997). 
The model is a long-wavelength adjustment to the FES94.1 hydrodynamic model 
of Le Provost et al (1994) (this is also the case for Richard Eanes's CSR3.0 
model). Cycles 9 to 71 of both TOPEX and POSEIDON data were used in the 
model development; an adjustment for the JGM-3 ephemerides was done so that 
the model is consistent with the second-generation orbits available for T/P. Model 
SR9601 04 differs from SR950308 only in some slight adjustments toward higher 
latitudes, which were necessary to extend the model into latitudes above 66°. 

The model is given on a 0.5° grid. 

The SR9601 04 ocean tide model has 7 primary tidal constituents: Oi, Pi, Ki, N 2 , 
M 2 , S 2 , and K 2 . The Qi constituent has been adopted directly from model 
FES95.2.1 of Le Provost et al. (1 994); this Qi is an assimilation solution, based 
on T/P data. In addition to these 8 constituents, 1 6 minor constituents are 
accounted for in the tidal-height predictions; these terms are inferred from the 
major tides by the admittance method (Munk and Cartwright, 1 966). Nodal 
correctipns for lunar tides are, of course, incorporated. Finally, the long-period 
tides are handled similarly to the T/P GDRs: they are assumed to be in 
equilibrium; the 15 largest spectral lines in the Cartwright and Edden (1973) 
tables have been used, including the 18.6-year nodal tide but not the permanent 
tide. 

In the polar oceans (latitudes above 66°), including Hudson Bay, the ocean-tide 
model is FES94.1 , although with some slight adjustments near the transition to 
SR960104 in order to allow for (relatively) smooth corrections without boundary 
jumps. Therefore, all latitudes overflown by the ERS-1 altimeter have ocean-tide 
corrections available. However, they should be considered less accurate than 
latitudes overflown by T/P where the SR9601 04 model directly incorporates the 
T/P measurements. 

The available global T/P-based tide model is not accurate in shallow seas. In 
these areas, the tides are spatially complex and the spatial averaging that is 
necessary to overcome tidal aliasing problems in altimetry cannot be performed 
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without considerable distortion of the tidal signal. The use of local hydrodynamic 
models may therefore often be preferable. 

The Pathfinder project has begun an effort to supplement the above global tide 
model with a number of local high-resolution hydrodynamic models for various 
marginal and inland seas. For the ocean-tide corrections in current use, the list 
below shows which model has been used in which location. 

• Mediterranean Sea : From Ganceil, Agelou, and Vincent (JGR, in press). 
Available on a 0.1° grid. 

• Persian Gulf : From Proctor, Flather and Elliott, Continental Shelf Research, 

1 4, 531 , 1 994. Available on a 5’ grid. 

• Gulf of Maine : From A. Lambert (Bull d'lnfor Mar T err, 1 1 0, 801 7, 1 991 ) 
which is a synthesis of results by Godin fCotidal Charts for Canada,' Dept. 
Fish. Oceans, 1 980) and Greenberg (Marine Geodesy, 2, 1 61 , 1 979). 
Includes the lower Bay of Fundy. Given on a 0.25° grid. 

• Gulf of St Lawrence : From A. Lambert (Bull d'lnfor Mar T err, 1 1 0, 801 7, 

1 991 ), essentially the same as Godin ('Cotidal Charts for Canada'). Given on 
a 0.25° grid. 

In addition, inland seas as well as the Black Sea and Baltic Sea are assumed to 
have no tides, which is not strictly true, but no accurate models are currently 
available. The one body of water with known large tides but with no available 
(and reasonably accurate) model is the Red Sea. 

We are accumulating additional local models for this project, and new versions of 
our altimeter products will gradually reflect this. 

The utility of local models for handling the tide corrections in marginal and shallow 
seas can be seen from the following experiments. Two-cycle differences of T/P 
collinear sea-surface heights were processed along three tracks segments(track 
1 01 across the Gulf of Maine and tracks 3 and 54 across the Persian Gulf). The 
rms of the height differences were tabulated using (a) no tide correction at all (b) 
the tide correction adopted for the second release of the GDRs (i.e. CSR3.0) and 
(c) the local models adopted in this project. 


• Table 4 Impact of local tide models: RMS of collinear differences (cm). 


Track 

No tidal correction 

GDR model 
correction 

Pathfinder local 
model correction 

101 

133 

95 

18.7 

3 

39 

34 

15.3 

54 

68 

75 

20.4 


The reduced variances resulting from using the local models provides convincing 
evidence of the value of these models. 
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Figure 13 Geographic distribution of the coverage of the NASA Ocean 
Pathfinder Project Ocean Tide Model. 
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Schrama-Ray 960104, a long-wavelength adjustment to FES94.1 
FES94.1 , a hydrodynamic model by Le Provost et al. [1994] 
Various local hydrodynamic models 



• Figure 14 Geographic distribution of the M2 ocean tide amplitude (cm) 
according to model SR950308. White lines denote cotidal phase lines. 
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Figure 15 Geographic distribution of the Pathfinder Ocean tide model 
differences with respect to CSR model (cm). 



Figure 15 shows the difference in cm between the Pathfinder ocean-tide model 
and the CSR3.0 model produced by Richard Eanes (University of Texas). The 
figure gives the maximum difference observed between the two model over a 
particular day (10 Oct 1995). (The band at the highest latitudes has been 
purposely increased— see below). These two T/P-based models are clearly very 
similarly over nearly over nearly all the deep ocean. In shallow and marginal 
seas, where the tides become larger and spatially complex, there are much larger 
differences between models. 

In the Pathfinder collinear data products, this image has been used for setting the 
“Suspect Tide Correction” flag bit. Wherever the tide-model differences exceed 5 
cm, this bit has been turned on. The latitudes outside T/P’s inclination (i.e. those 
exceeding 66 degrees) require special handling. In those latitudes, both the 
Pathfinder model and CSR3.0 use the hydrodynamic model (FES94.1) of Le 
Provost, which is, of course, much less reliable than the T/P-based model. To 
reflect this, the model differences in these regions have been artificially inflated, 
and the suspect tide bit is always on. 

2.5.3. Load Tide Correction 

The load tide describes the local vertical displacement of the solid earth caused 
by the weight of the overlying ocean tide. The load tide must therefore be 
consistent with the mode! used for the ocean tides. The load tide is here modeled 
as a purely elastic response to ocean loading, using a high-degree expansion in 
spherical harmonics (to degree and order 360). The response of each degree 
depends on the loading Love number h'„ ; these have been adopted from the 
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calculations by Farrell (1972). Further detailed descriptions of this method of 
computing the load tide can be found in Ray and Sanchez (1989). The amplitude 
of the M 2 load tide, in millimeters, is shown in Figure 1 6. The largest amplitudes 
reach 5 cm off the coast of Brazil. For the diurnal load tides, the largest 
amplitudes are in the northern Pacific (e.g. the Gulf of Alaska) and off the coast 
of Antarctica; in these regions the K-i load tide reaches nearly 3 cm. 

• Figure 16 Geographic distribution of the M 2 amplitude of the ocean load 
model (mm). 



2.5 A Pole Tide Correction Version 2.0 

The pole tide is the response of the ocean (and also of the elastic earth) to 
variations in the centrifugal force caused by wobbling of the earth's rotation axis. 
The pole tide has two dominant frequencies: annual and 14-month, the latter 
being the period of the Chandler Wobble. 

The pole tide correction to the altimetry is a geocentric correction, meaning that it 
includes both the ocean and the solid-earth pole tide. An equilibrium response 
has been assumed for the ocean. The computed effect is given by Munk and 
MacDonald, (1 960): 
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Y (mas) X (mas) 


Eq. 6 Pole Tide elevation = A *sin(2*lat) * [x*cos(lon) + y*sin(lon)] 


Where, x and y denote the position of the pole along the Greenwich and 90°E 
meridians, respectively. Because the pole position has a non-zero mean, which 
itself drifts slowly with time, a linear trend is first removed from (x, y). The (x, y) 
positions of the pole are obtained from the standard 5-day time series of the 
International Earth Rotation Service, Paris. 

The amplitude A = -68.85 mm when (x, y) are expressed in arcsec. However, 
when the altimeter is overflying an inland sea, not connected to the global ocean, 
only the solid-earth pole tide is effective. For that case, A is scaled by the Love 
number factor h/(1 + k) = 0.47. 

♦ Figure 17 The (x, y) data of the Pole tide model since 1958 ; including the 
linear trend. 
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Latitude 


Figure 18 The geographic distribution of the pole tide correction (mm) for 
ERS-1 Phase C cycle 14. 
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2.5.5. Crosstrack Geoid Correction 

Altimeter satellites are maintained in a repeating orbit to facilitate the separation of 
sea-height variations from the geoid. Perturbations in the satellite orbit result in 
excursions from the "exact repeat" groundtrack. For Geosat this drift caused the 
groundtrack to vary by more than ±1 km about the nominal repeat path. This 
misalignment leads to an error in the estimates of sea surface height variations 
because of the local slope in the geoid (Brenner et al., 1 990). High resolution 
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mean sea surfaces have been constructed from the combined Geosat, ERS-1 , 
and T/P altimeter data sets enabling the computation of accurate geoid slopes 
norma) to the reference repeat orbit. The geoid gradients used in these 
PATHFINDER products were computed from the mean sea surface 
CSRMSS95IB, developed at the University of Texas. 

Does this correction demonstrably improve the data? An examination of the 
variances of sea-surface height differences after applying or not applying this 
correction (in the manner of Ray, Koblinsky, and Beckley, 1991) suggests that the 
correction is indeed of value. Using data from both Geosat and T/P, we found 
variance reductions for both satellites of order 1 cm 2 . 

• Figure 19 The geographic distribution of the ascending and descending orbit 
cross track geoid gradient correction (mm) for Geosat, ERM Cycle 27. 



2.5.6. Atmospheric Loading (Inverted Barometer) 

Some of the Pathfinder datasets have been adjusted for atmospheric loading. 
This is useful in many applications when the effects of atmospheric loading 
are not of interest, for example in the estimation of surface geostrophic 
currents. This loading ‘correction’ is done by use of a simple, isostatic, 
inverted barometer assumption (Gill, 1982): 
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• Eq. 7 Inverted Barometer Correction = - 9.948 * ( P - 1013.3) millimeters 

Where, P is the atmospheric surface pressure in millibars and 1013.3 is the 
global mean atmospheric pressure. In our data processing, these surface 
pressures are consistent with the ones used for the dry-troposphere 
correction. For T/P, Geosat, and ERS-1 the pressures used are the 
operational products of the European Centre for Medium-range Weather 
Forecasting (ECMWF). For Seasat, we have used surface pressures from 
the Fleet Numerical Ocean Center atmospheric model. 

• Figure 20 The geographic distribution of the atmospheric load (mm) 
determined from the ECMWF surface pressure fields and the inverse 
barometer response model. This is the mean load during cycle 2 of the ERS-1 
Phase C Exact Repeat Mission , i.e. during the 35-day interval starting 19 May 
1992. 
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3. Product 1: The Coliinear Data Set 


3.1. Description 

The Ocean Altimeter Pathfinder Project began its work by reproducing all 
GDRs for each mission with revisions based upon the corrections described 
above. These new GDRs formed a basis for the formation of ‘easy-to-use’ 
products. The first product created from these data and distributed to the 
community was the simplified coliinear data set. This data set consists of files 
containing regularly spaced, spatially-indexed, coliinear sea surface heights 
with respect to a reference mean sea surface. These data were derived from 
coliinear mission data only. For the temporal coverage of the coliinear 
missions see Table 5. 


• Tables Coliinear data coverage. 


YEAR 

Jan 

Feb 

Mar 

El 

May 

June 

July 

Aug 

Sep 

Oct 

Nov 

Dec 

78 







2 

2 

2 
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Geosat 

Geosat 

87 

Geosat 

Geosat 

Geosat 

Geosat 

Geosat 

Geosat 

Geosat 

Geosat 

Geosat 

Geosat 

Geosat 

Geosat 

88 

Geos at 

Geosat 

Geosat 

Geosat 

Geosat 

Geosat 

Geosat 

Geosat 

Geosat 

Geosat 

Geosat 

Geosat 

89 

i 


1 

i 

1 

i 

1 

1 1 

1 — 

1 

ESI 

i 

90 
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"" 3 

3 

3 

3 

" ~ ~3 " 

92 

Bi 

Bi 

3 

3 

ERS-1 

ERS-1 

ERS-1 

ERS-1 

ERS-1 & 
T/P 

ERS-1 & 
T/P 

ERS-1 & 
T/P 

ERS-1 & 
T/P 

93 

ERS-1 & 
T/P 

ERS-1 & 
T/P 

ERS-1 & 
T/P 

ERS-1 & 
TP 

ERS-1 & 
T/P 

ERS-1 & 
T/P 

ERS-1 & 
T/P 

ERS-1 & 
T/P 

ERS-1 & 
T/P 

ERS-1 & 
T/P 

ERS-1 & 
T/P 

ERS-1 & 
T/P 

94 

mrm 

Ball 

WSSMi 

TP 

T/P 

T/P 

T/P 

T/P 

T/P 

T/P 

T/P 

T/P 

95 

TP 

T/P 

TP 

ERS-1 ,2 
& TP 

ERS-1, 2 
& TP 

ERS-1 ,2 
& T/P 

ERS-1 ,2 
& T/P 

ERS-1, 2 
& T/P 

ERS-1, 2 
& T/P 

ERS-1, 2 
& T/P 

ERS-1, 2 
& T/P 

ERS-1, 2 
& T/P 

96 

ERS-1 ,2 
& TP 

ERS-1. 2 
& TP 

ERS-1 .2 
& T/P 

ERS-1, 2 
&T/P 

ERS-1, 2 
& T/P 

ERS-2& 

T/P 

ERS-2 & 
T/P 

ERS-2& 

T/P 

ERS-2& 

T/P 

ERS-2& 

T/P 

ERS-2 & 
T/P 

ERS-2 & 
T/P 

97 

ERS-2& 

TP 

— — — r 1 

ERS-2& 

TP 

--GEOSAT C 

ERS-2& 

T/P 

Esa In i ytiy n 

ERS-2& 

T/P 

t WKssecr 

ERS-2 & 
T/P 

Kcauss u! a 

ERS-2& 

T/P 

Missive orbit < 

ERS-2& 

T/P 

TOT 

ERS-2& 

T/P 

ERS-2 & 
T/P 

ERS-2 & 
T/P 

ERS-2 & 
T/P 

ERS-2 & 
T/P 


2 - Seasat never in exact repeat orbit 

3 - ERS-1 3-day repeat data not processed for pathfinder 


Altimeter satellites that are maintained in a repeating orbit facilitate the 
separation of sea height variations from the geoid. The term ‘collineari refers 
to sea surface heights for a particular ‘exact repeat orbit’ mission that have 
been georeferenced to a specific groundtrack. The coliinear data file contains 
sea surface heights for each orbit cycle at fixed locations thus allowing for the 
direct analysis of Sea Surface Height (SSH) variability. 
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The sea surface heights were corrected for all geophysical, media, and 
instrument effects are given at 1 second intervals (about 7 km, less for T/P) 
along the reference track using the algorithms described in Section 2. For 
each satellite files were created that contain: 

a) georeference locations; 

b) reference mean sea surface height relative to the reference ellipsoid 
described in Table 3 at each georeferenced location; 

c) the time series of corrected sea surface height at each location; 

d) a quality control file that contains a flagword for each observation and 

e) a file that provides the observation time at index 1 for each revolution. 
Software to access all of the files is provided by the project (see below). 

3.2, Processing Procedure for Coilinear data product 

• Obtain, analyze, and apply the latest state-of-the-art corrections and models, 
consistent across missions, whenever possible. 

• Compute improved GDRs. 

• Spatially and temporally interpolate to reference ground-track at 1/sec 
sampling 

• Apply editing criteria to extract bad data 

• Compute flag data based on quality control criteria 

3.3. General Format for Sea Surface Height ( SSH ) File 

The format of the sea surface height file is a sequence of binary records that 
contain: 

Alongtrack Index 

Revolution Number 

Array of Sea Surface Heights in mm 


All the variables are IEEE 2-byte integers. The array length depends on the 
number of repeat cycles for the particular database version. 



Table 6 Summary Information of Altimeter Repeat Missions for Ocean 
Pathfinder Data Sets as of 12/1/97. 


Satellite 



Revolution 
Period in 
Seconds 


Number of 
Cycles 


Span of Data 


TOPEX/ 

POSEIDON 

9.92 

127 

6745 

154 

92/09/23 - 
96/11/28 

ERS-1 
Phase C 

35.00 

501 

6038 

18 

92/04/14 • 
93/12/20 

Geosat 

(ERM) 

17.05 

244 

6038 

62 

86/11/08- 

89/09/30 


3A Dataset File Names and Descriptions 

ssh.dat File containing the sea surface height residuals (mm) with along 

track index and revolution number 

ssh_noib.dat File containing the sea surface height residuals (mm) with along 
track index and revolution number. Sea surface heights DO 
NOT have the inverted barometer correction applied 

reforb.dat File containing geodetic North latitude and East longitude in 
microdegrees of georeferenced locations 

mss.dat CSR mean sea surface (CSRMSS95IB) interpolated to the 
indexed georeferenced locations (cm) 

flagword.dat Data quality flagword for each observation based on certain 
criteria. 

time.dat Observation time at index 1 for each revolution. All times are in 
terms of Modified Julian Dates. 

directry.dat Directory of direct-access record numbers referenced by 
revolution and alongtrack index . 

Please note that a revolution is defined as the ground track between consecutive 
ascending equator crossings. The SSH array contains sea surface heights 
referenced to the University of Texas Center for Space Research (CSR) mean 
sea surface (CSRMSS95IB). The mean sea surface was computed from the 
combined Geosat, ERS1 and T/P altimeter data sets. The Geosat and ERS1 
data were globally and simultaneously crossover-adjusted with respect to 
TOPEX/POSEIDON. 

Note: It was discovered in 1 996 that a miscalculated instrument correction was 
applied to the TOPEX range measurement resulting in an ~13 cm bias. 
Therefore, the mean sea surface is elevated 1 3 cm more than it should be with 
respect to the reference ellipsoid. 
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3.5. Calculated Offsets between Missions 


Mission offsets were computed by solving for Ar, Ax, Ay, and Az from crossover 
residuals with respect to a T/P 1993 mean annual profile. Annual mean profiles 
were computed for Geosat and ERS-1 and crossed into T/P. 
i.e. : 

1993 ERS-1 mean profile x T/P 1993 mean profile 
1 987 Geosat mean profile x T/P 1 993 mean profile 

Table 7 Offsets (cm) from crossover residuals relative to TOPEX/POSEIDON . 


Satellite 

Ar 

Ax 

Ay 

AZ 

Geosat 

12.4 

7.9 

3.3 

1.8 

ERS-1 

-43.5 

-2.5 

-2.5 

-1.5 


In order to correct the collinear data for these offsets the following formula should 
be used: 

• Eq. 8 Corrected SSH = SSH + (Ar+Ax*cos(A)*cos(<j))+ Ay*cos(X)*sin((f>)+ Az*sin(X)) 
Where, $ = East longitude, North latitude 

The mean sea surface is provided so that sea surface height can be referenced 
to the TOPEX standard reference ellipsoid, if desired, see Table 3. 

3.6. Updates to Collinear Data Set 

3.6.1. TOPEX/POSEIDON 

TOPEX/POSEIDON Version 2.0 Updates from Version 1.0: 

Revised Feb. 19, 1997: Version 2.0 

NASA/GSFC Pathfinder TOPEX/POSEIDON Georeferenced Altimetry 
Heritage: PODAAC GDR5.0 for TOPEX / PODAAC MG DR for POSEIDON 

A secular trend for the Earth's mean pole position has been removed before 
computing the pole tide. Data over inland seas have been restored. (The 
exception is the Red Sea, for which reliable tide corrections are lacking.) Ocean 
tide quality flag bit (see flag documentation). Sea surface height file contains 
cycles 1-154, with the exception of Poseidon cycles 138 and 150. 

Sea surface heights were computed with JGM3 orbits. The precise orbital 
ephemerides for T/P are identical to those used in the second version of the T/P 
GDRs, as computed by the Space Geodesy Branch at NASA/GSFC. These 
ephemerides are computed for nominal 10-day arc segments, corresponding to 
the groundtrack repeat period of the satellite. They are based on force models 
derived from the JGM-3 earth gravity model, a T/P-based ocean-tide model, a 
sophisticated 'box-wing' satellite model for drag and radiation forces, and on laser 
and doppler tracking data. Nerem et al. (1993) and Tapley et al. (1 994a) describe 
earlier versions of the orbits. Those papers give many further details on the 
ephemeris calculations. The present version of the orbits is described by 
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Marshall et al. (1995), who examine in detail the error budget while presenting 
strong evidence for a radial orbit accuracy of about 3 cm rms. 

Each record contains an array of sea surface heights for the 154 ten day repeat 
cycles at a single location defined by a revolution number and along track index: 

* Table 8 TOPEX/POSEIDON Parameter File Record Description. 


Variable 

Format 

Word Size in Bytes 

Along Track Index 

1*2 

2 

Revolution Number 

1*2 

2 

Sea Surface Residual 
Height Array (mm) 

1*2 

308 


Editing was performed on the data based on these criteria. A flagword for each 
sea surface residual height was set based on certain tests. 

3.6.2. ERS-1 Version 2.0 Updates front Version 1 .0: 

Revised Feb. 19, 1997: Version 2.0 

NASA/GSFC Pathfinder ERS-1 Phase C Georeferenced Altimetry Fleritage: 
CERSAT OPR 3.0 

A secular trend for the Earth's mean pole position has been removed before 
computing the pole tide. Data over inland seas have been restored. (The 
exception is the Red Sea, for which reliable tide corrections are lacking.) Major 
update: Sea surface height is based on Delft DGM-E04 orbits with 1 .7 ms time 
tag bias applied, (see Section 2). Ocean tide flag bit (see flag documentation). 

• Table 9 ERS-1 Parameter File Record Description. 


Variable 

Format 

Word Size in Bytes 

Along Track Index 

1*2 

2 

Revolution Number 

1*2 

2 

Sea Surface Residual 
Height Array (mm) 

1*2 

36 


Editing was performed on the data based on these criteria. A flagword for each 
sea surface residual height was set based on certain tests. 

3.6.3. Geosat Version 2.0 Updates from Version 1 .0: 

Revised Feb. 19, 1997: Version 2.0 

NASA/GSFC Pathfinder Geosat ERM ( Exact Repeat Mission ) Georeferenced 
Altimetry 

Fleritage: NOAA/NODC GEMT2-Based GDR (with many improvements, 
including JGM-3 orbits). 
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A secular trend for the Earth's mean pole position has been removed before 
computing the pole tide. Data over inland seas have been restored. (The 
exception is the Red Sea, for which reliable tide corrections are lacking.) Wet 
tropospheric correction based on NCEP reanalysis model (PWAT). (AMS 
Bull.77#3,3/96.) (see Section 2 ). Sea state bias correction. (Gaspar, 1996) ( see 
Section 2 ) Ocean tide quality flag bit ( see flag documentation ) 


• Table 10 Geosat Parameter File Record Description. 


Variable 

Format 

Word Size in Bytes 

Along Track Index 

1*2 

2 

Revolution Number 

1*2 

2 

Sea Surface Residual 
Height Array (mm) 

1*2 

124 


Editing was performed on the data based on these criteria. A flagword for each 
sea surface residual height was set based on certain tests. 

3.7. Validation 


The validation of the colinear products are summarized in Table 11 (RMS 
companion). This validation is fully described in Report #2: Validation Handbook. 


• Table 1 1 Global Statistics of Sea Surface Height Variations Computed from 
Collinear Altimetry Compared to Tide Gauge Measurements. 



# of Stations 

Mean RMS (cm) 

Mean Correlation 

T/P 

76 

6.77 

0.71 

ERS-1 Phase C 

78 

6.47 

0.70 

Geosat ERM 

58 

8.55 

0.62 




























4.1. Description 

A time series of gridded monthly sea surface height variations with respect to a 
1 993 mean reference has been generated from all available altimeter 
observations provided by Seasat, Geosat, ERS-1 , and T/P. The grids are 
created from height anomalies of cross over or collinear data relative to a mean 
reference set of collinear ground tracks. The time period for these grids is 
described in Table 12. 


Table 12 Summary of Monthly Coverage of Grids between 1978 and 1997. 


Year Jan 


78 


79 


Feb I Mar 


y 

Au 

g 

Se 


Oct Nov ! Dec 




Geosat 

Geosat 

Geosat 

Geosat 

Geosat 

Geosat 

Geosat 

Geosat 

Geosat 

Geosat 

Geosat 

Geosat 

Geosat 

Geosat 

Geosat 

Geosat 

Geosat 

Geosat 

Geosat 

Geosat 

Geosat 

Geosat 

Geosat 

Geosat 

Geosat 

Geosat 

Geosat 

Geosat 

Geosat 

Geosat 

Geosat 

Geosat 

Geosat 

Geosat 

Geosat 

Geosat 



92 


~=-S ' 

m 

|| 

ERS-1 

ERS-1 

ERS-1 

ERS-1 

ERS-1 

ERS-1 & 
T/P 

ERS-1 8c 
T/P 

ERS-1 8c 
T/P 

93 

ERS-1 & 
T/P 

ERS-1 & 
T/P 

ERS-1 & 
T/P 

ERS-1 & 
T/P 

ERS-1 & 
T/P 

ERS-1 & 
T/P 

ERS-1 & 
T/P 

ERS-1 & 
T/P 

ERS-1 8c 
T/P 

ERS-1 & 
T/P 

ERS-1 & 
T/P 

ERS-1 8c 
T/P 

94 

ERS-1 & 
T/P 

ERS-1 & 
T/P 

ERS-1 & 
T/P 

ERS-1 & 
T/P 

ERS-1 & 
T/P 

ERS-1 & 
T/P 

ERS-1 & 
T/P 

ERS-1 8c 
T/P 

ERS-1 8c 
T/P 

ERS-1 8c 
T/P 

ERS-1 8c 
T/P 

ERS-1 8c 
T/P 

95 

ERS-1 & 
T/P 

ERS-1 & 
T/P 

ERS-1 & 
T/P 

ERS-1 ,2 
&T/P 

ERS-1 ,2 
8c T/P 

ERS-1, 2 
& T/P 

ERS-1 ,2 
& T/P 

ERS-1, 2 
8c T/P 

ERS-1, 2 
8c T/P 

ERS-1 ,2 
8c T/P 

ERS-1, 2 
8c T/P 

ERS-1, 2 
8c T/P 

96 

ERS-1, 2 
& T/P 

ERS-1, 2 
& T/P 

ERS-1, 2 
& T/P 

ERS-1 ,2 
&T/P 

ERS-1 ,2 
& T/P 

ERS-2 8c 
T/P 

ERS-2 8c 
T/P 

ERS-2 8c 
T/P 

ERS-2 8< 
T/P 

ERS-2 8c 
T/P 

j ERS-2 8c 
T/P 

ERS-2 8c 
T/P 

97 

ERS-2 8c 
T/P 

ERS-2 8c 
T/P 

ERS-2& 

LJE* 

ERS-2 8c 
T/P 

ERS-2& 

T/P 

ERS-2& 

T/P 

ERS-2 & 
T/P 

ERS-2 8c 
T/P 

ERS-2 & 
T/P 

ERS-2 8c 
T/P 

ERS-2 8c 
T/P 

ERS-2 8c 

. Iff 



4*2. Processing Procedure for Grid Data Product 

• Obtain, analyze, and apply the latest state-of-the-art corrections and models, 
consistent across missions, whenever possible. 

• Compute improved G DRs. 
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• Spatially and temporally interpolate to reference ground-track at 1/sec 
sampling. 

• Apply editing criteria to extract bad data. 

• Compute mean reference surface from mean T/P and ERS data sets. 

• Compute crossover anomalies for non-ERS or T/P data sets. 

• Adjust cross over or collinear anomalies into mean su rface. 

• Compute grid from anomalies within a given month. 

4*3. Algorithm 

4*3.1. Mean Reference Surface 

The mean reference surface was computed from merged ERS-1 and T/P annual 
mean heights along ground-tracks for 1993 by adjusting the 35-day ERS-1 mean 
profile for 1993 into the T/P geodetic frame. The ERS-1 Phase C altimetry was 
first updated with Delft DGM-E04 JGM3 based orbits and a number of improved 
geophysical and altimeter range corrections (see Section 2). A mean profile is 
then computed at each geo-referenced location employing standard edit 
procedures, including a minimum cycle requirement to avoid seasonal aliasing. 
Adjustment of the ERS-1 mean profile with respect to the T/P mean profile was 
performed by crossover minimization within 30° zonal bands modeled by a simple 
tilt and bias. Crossover analysis of the adjusted ERS-1 mean profile with the T/P 
mean resulted in a residual rms of less than 3 cm. The " blending " of ERS-1 and 
T/P mean profiles generates a dense ground-track providing a statistical robust 
sample of crossover and collinear sea surface height residuals (see Figure 21). 
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• Figure 21 Coverage of blended ERS-1 and TIP groundtracks that form 
reference surface for grids. 



T/P- 


ERS— 1 


4.3.2. Adjustment Procedure 

Significant improvements have been made to the orbit computation based on the 
JGM3 gravity field for Seasat, Geosat, and ERS missions. However, the orbits 
have not reached the accuracy of the T/P altimeter mission. Therefore, sea 
surface height variations computed from Seasat, Geosat, and ERS crossover or 
collinear residuals with respect to the 1993 reference have been adjusted to 
remove systematic orbit errors and geocentric offsets. These errors are modeled 
with a 1 cycle per revolution sinusoid and bias term from the simultaneous 
crossover minimization of single and dual crossovers to the T/P mean profile. 

Reduction of systematic orbit errors employing crossover minimization techniques 
has a long heritage (e.g., Rapp, 1977; Tai, 1988, 1994; Hwang, 1989; Yi, 1995). It 
is well known that the orbit error has dominant power at one cycle per revolution 
(cpr) with secondary peaks at 2 cpr and 3 cpr (Rosborough and Tapley, 1 987; 
Engelis, 1988; Chelton and Schlax, 1993). A low degree polynomial or Fourier 
series have been used to model the orbit error. An important aspect in using 
these kinds of error models is that improper use of the error model could remove 
ocean variability (Tai, 1 989; Wagner and Tai, 1 994). In order to minimize the loss 
of ocean dynamics, the polynomial expansion or Fourier series of the lowest 
possible degree and order have to be chosen. In our study the sinusoidal error 
model is chosen for one revolution of the altimeter data: 

(1) r(t) = a +b cos cof + c sin cof 


45 


where o)=2k/T is the one cpr for the nodal period of the altimeter mission, a, 
b and c and are the estimated coefficients for a single revolution. This error 
model has an along-track wavelength of 40,000 km, thus any effect on the ocean 
variability has to be at signals with wavelengths of 40,000 km or longer. 

The residual sea surface at the single-crossover point (single altimeter mission) 
is: 

(2) Ah, r r(ti)-r(ti)+ Vfatj) 

where t„tj are the times of the ascending and descending arcs at the crossover 
point; 5 (&0 is the time varying sea surface height and error of the environmental 
corrections at the crossover point. Assuming the TOPEX/POSEIDON mean 
profile is free of orbit error, the residual sea surface height at the dual crossover 
point (dual altimeter missions) is: 

(3) A/?/y= r(tj) + 8* 


Where 5 is the sea surface height variation at the crossover point respect to 
the TOPEX/POSEIDON mean profile and the error of the environmental 
corrections in both altimeter missions. 

Equation (2) and (3) form two observation equations, and the parameters a, b 
and cfor each revolution are estimated through the standard least squares 
adjustment. 

The crossover adjustment procedure has been applied to the Seasat, ERS-1 
phase C and Geosat missions. The following table gives the RMS values of the 
residual sea surface height at the crossover point before and after adjustment. 
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Table 13 RMS values of the residual sea surface height at the crossover 
points before and after the adjustment (meters). 



Single Satellite 

E 

Dual Satellites 

Number 



■ 

Number 

mmsm 

IfcaiKIlUi 

Seasat 

(3days) 

Cycle 2 

516 

0.265 

0.137 

1 

3797 

0.534 

0.136 


EEEIM 

521 

0.257 

0.127 


3728 

0.509 

0.128 


ksssieh 

461 

0.329 

0.121 


3556 

0.578 

0.139 










Seasat 
(17 day) 

Cycle 1 

3050 

0.523 

0.121 

. 

1 

9547 

0.558 

0.134 



6562 

0.312 

0.119 


14063 

0.536 

0.135 










Geosat 

ERM 

Cycle 21 

14231 

0.151 

0.101 

1 

20862 

0.197 

0.120 



17218 

0.163 

0.088 


25066 

0.212 

0.105 



9736 

0.222 

0.101 


17029 

0.226 

0.110 










Geosat 

GM 

Cycle 14 

35976 

0.153 

0.086 

1 

37732 

6.181 

0.098 



22255 

0.136 

0.084 


30302 

6.181 

0.100 


Cycle 25 

19081 

0.128 

0.084 


27196 

0.171 

0.102 










m&Emm 

Cycle 10 

32507 

0.126 

0.093 

1 

49441 

0.421 

0.088 



28467 

0.129 

0.097 


35839 

0.475 

0.100 



27109 

0.117 

0.093 


39075 

0.449 

0.099 


4.3.3. Grid Computation 

Monthly grids were derived by estimating at each point P on a 1° X 1° regular grid 
a weighted average of adjacent (within a 3° radius) cross over or collinear 
differences with respect to the mean 1 993 reference surface. Figure 22 shows 
typical crossover density of points between Geosat GM and the reference surface 
in mid-latitude regions. In lieu of using an optimal estimator based on local values 
of the height covariance for smoothing weights, the data weights were derived 
using an isotropic Gaussian function depending only on the distance from the 
calculation point P(<j>A) to the data point at (fah). 
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• Figure 22 Groundtrack coverage with respect to 1993 reference frame from 
a) Geosat Primary Mission groundtracks and b) Geosat Exact Repeat Mission 
ground tracks. 



T/P ERS-1 GEOSAT GM 



T/P ERS-1 GEOSAT ERM 
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Latitude 


• Eq. 9 Grid Algorithm 


h P 0 , A ) = 


I A,) 


i 


W 


Where, wp exp(-a d 2 ); d - spherical distance between P(0,A.) and the data points 
((j)i,Xi); a - smoothness parameter defined by the half weight x, such that: 
exp(-crr 2 ) =1/2 a=ln(2)x' 2 . For our grids, we have selected x=1 °. Grids of the 

sea surface height variations are computed using an interpolation procedure to 
predict a weighted average SSH residual within a defined search radius of the 
desired grid node. The procedure is fully described in Nerem et al. (1 994). . 
Several tests showed that the crossover densities were sufficient to permit a half 
weight parameter of 1°, and a search radius of 3° (see Figures 23 & 24). The 
validation of these grids against tide gauge height changes (summarized in Fig. 
27) showed that these parameters represent the original data on a global scale. 
Report #2: Validation Handbook provides complete validation details regarding 
these grids. 

• Figure 23 Global density of crossover measurements within each grid node 
from Geosat Primary Mission fora month (November 1985). 
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Latitude 


Figure 24 Gridded crossover sea surface height residual for Geosat GM 
November 1985 with respect to 1993 reference. 
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Figure 25 Validation of Gridded Altimetry Against Tide Gauge Network . The 
Mean RMS Difference of the Monthly Grids vs. Tide Gauge is 4.6cm. 
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